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Context. For Galactic open clusters, fundamental parameters such as age or reddening are usually determined independently of their 
integrated colours. For extragalactic clusters, on the other hand, they are derived by comparing their integrated colours with predictions 
of simple stellar population (SSP) models. 

Aims. We search for an explanation of the disagreement between the observed integrated colours of 650 local Galactic clusters and 
the theoretical colours of present-day SSP models. 

Methods. We check the hypothesis that the systematic offsets between observed and theoretical colours, which are (B—V) ~ 0.3 
and (J—K s ) ~ 0.8, are caused by neglecting the discrete nature of the underlying mass function. Using Monte Carlo simulations, 
we construct artificial clusters of coeval stars taken from a mass distribution defined by an Salpeter initial mass function (IMF) and 
compare them with corresponding "continuous-IMF" SSP models. 

Results. If the discreteness of the IMF is taken into account, the model fits the observations perfectly and is able to explain naturally 
a number of red "outliers" observed in the empirical colour-age relation. We find that the systematic offset between the continuous- 
and discrete-IMF colours reaches its maximum of about 0.5 in (B—V) for a cluster mass M c = 10 2 mg at ages logt ~ 7, and 
diminishes substantially but not completely to about one hundredth of a magnitude at logt > 7.9 at cluster masses M c > 10 5 rriQ. 
At younger ages, it is still present even in massive clusters, and for M c ^ 10 4 m© it is larger than 0.1 mag in (B—V). Only for very 
massive clusters (M c > 10 6 mg) with ages logt < 7.5 is the offset small (of the order of 0.04 mag) and smaller than the typical 
observational error of colours of extragalactic clusters. 
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• 1 . Introduction respect to the observations for clusters younger than log t < 8.5. 

O . .... This is both remarkable and important, since most data such as 

Using data on accurate and homogeneous spatio-kinematic- and m collected nowad for the bulk of extragalactic 

^ , photometric membership for 650 Galactic open clusters we duste are derived b comparing their i ntegrated ligh t with 

g ■ previously compute d their integrated magnrtu des in B, V, J,H, ^ ed ictions of SSP models (see Bik et al 2003) The 

^ , and A s -passbands ( |Kharchenko et al. | |2^ The magnitudes modelfj m duced b evolutionary synthesis codes simulating 

> . are based on accurate and uniform data from the ASCC-2.5 cat- the exte nsive and complicated stellar populations of galaxies. 

£j ; alogue and were computed by adding the individual luminosities However> th ^ also believed t0 be a p propriate for describing 

, of the most secure cluster members. In co ntrast to previous lists dugter lations more than five orders of magnitude less 

^ ■ of integra t ed m agnitudes (e.g., those of | Battinelli et al. || l994t massiye (and fl far si le r). Although star clusters 

| Lataetal. | |2002j of Galactic star clusters, our data provide an fee influenced b effects that ^ irre levant to galaxies 

independent and more importantly uniform dataset. Hence, our ( bw number because of the discreteness of the 

integrated magnitudes can and should be used as benchmarks for feal ^ Qr d amical mass . loss by evaporation of stars), the 

studies of the populations of extragalactic star clusters. yalidi of ^ ssp ch for staf dusters cafl bg tested from 

Based on these uniform data we found disagreement ge neral considerations ( |Cervifio & LuridianalllOOl . In contrast, 

between the observed colours of our cluster sample and b ^ wg j ^ issue of ^ discreten e SS G f the IMF 

theoretical colours derived from simple stellar population to the spe cific case of Galactic open clusters. 
(SSP) models. The comparison of the Kharchenko et al. 



2009 j) data with present-day S S P mod els such as GALEV A_colour_offset of this kind was found previously by 

Anders & Fritze-v. Alvenslebenl l200l and Starburst99 |Laia^t_aL| (I200J), in the colour indices U-B, B-V, V-R, and 

Vazquez & Leithererl l2005h Illustrated a substantial offset of V-I. However, they accept that only the offsets in V-R and 

the order of 0.2 - 0.3 mag in (B—V) of the model colours with V-I represent significant deviation of their observations from 

the model predictions. We suspect that the discreteness of the 

Send offprint requests to: R.-D. Scholz IMF in open clusters may explain these discrepancies. To study 
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Figure 1. Observed colours of 650 Galactic open clusters com- 
pared to theoretical colours computed from standard SSP mod- 
els. The upper row is for B—V, the bottom one for J—K s . The 
left column shows "integrated colour vs. log t" diagrams, the 
right one the distributions of colour indices. Dots show Galactic 
open clusters from our sample. Black dots are young clusters 
used for constructing the colour distribution, grey dots the older 
clusters. The curves are the SSP model tracks (the solid one is 
from Starburst99, the dashed one is from GALEV).The observed 
colour distributions are shown as the filled cyan histogram, the 
GALEV model distribution as the hatched red one. 



this issue in a more systematic way, we constructed our own SSP 
model that takes into account the effect of low number statistics 
(discreteness), which is a common problem for open star clus- 
ters. We cross-checked a continuous version of our code with 
GALEV and found good agreement. 

The discreteness of the IMF itself is a natural assumption for 
any stellar ensemble consisting of individual objects. However, 
when one considers vast stellar populations (e.g., galaxies) that 
densely populate the entire colour-magnitude diagram, one can 
adopt a continuous IMF, which is more convenient for different 
technical reasons. For small populations (stellar clusters), the as- 
sumption of a discrete IMF is clearly appropriate.. This letter de- 
scribes the solution of a restricted problem, namely the influence 
of the discrete-IMF approach on the colours of young clusters. 
The full scope of results related to the IMF-discreteness effect 
on star clusters will be reported elsewhere. 



2. The SSP models and observations 

For the determination of integrated magnitud es (and colours) of 
Galact ic open clusters, we refer the reader to iKharchenko et all 
(2009), which also gives an overview of the cluster sample and 
references for further reading. 

For simplicity, we refer to our SSP model hereafter as the 
"realistic" model, and call the literature models "standard" mod- 
els. The standard SSP models were computed with help of the 
online servers provided by the GALE\fl and Starburst9S0 sites. 



http : / /www . galev . org/ 
2 http : / /www . st sci . edu/ science/ starburst 99/ 



The input parameters of the realistic model are the total clus- 
ter mass M c , the initial mass function (IMF) of cluster stars 

/(to) = dN/dm , 

(N being the number of stars with masses between to, m+drri), 
and the lower and upper limits (to; and m u ) of the range of 
the allowed masses of cluster stars. We adopt the Salpeter IMF, 
/(to) = km~ a , a = 2.35, where k is a normalisation factor 
that depends on the assumed cluster mass. The upper mass limit 
is confined to the high-mass end of the adopted grid of evolu- 
tionary tracks/isochrones. Only at the very early stages of cluster 
evolution does m u have an impact on integrated colours. Given 
a fixed IMF slope, the lower limit governs the population of the 
massive end of the IMF at given M c . We consider M c and mi 
here as free parameters of the adopted model. The model stel- 
lar sample was generated from the IMF with a random number 
generator. 

The photometric properties of the model are defined by the 
implemented grid of the isochrones. For our model, we have 
used the grid provided by the Padova group via the online server 
CME0 According to our experience, the Padova isochrones fit 
the brighter parts of the open cluster CMDs reasonably well over 
a wide range of cluster ages. These parts are responsible for the 
integrated colours. We used a solar metallicity grid (Z = 0.019), 
and retrieved the passbands U, B, V, R, I and J, H, K s with an 
age range logt = 6.0 — 10.2 and a step in logt of 0.05. For 
compatibility, we selected the same set of parameters (i.e., un- 
derlying IMF, mass range limits, isochrone set) for the standard 
models. However, since the Padova group regularly updates their 
isochrones, full compatibility between the isochrone grids used 
by the realistic and standard models cannot be guaranteed. 

In Fig.Q] we compare the location of the observations with 
the curves of the standard model in the colour-age diagram. 
GALEV and Starburst99 tracks agree reasonably well with each 
other. The above-mentioned offset of the models with respect to 
observed colours for log t < 8.5 is obvious in both of the colour- 
age diagrams. 

The observed clusters show a clear pattern, which at log t < 
8.5 is represented by "main sequence" clusters (not containing 
red giants or supergiants). The considerable redward spread of 
observed points cannot be attributed to photometric errors; it is 
provided in this part of the diagram by a few clusters containing 
yellow or red supergiants. For logt > 8.5, one finds, in gen- 
eral, clusters with red giants. The colour spread of older clusters 
is caused by different evolutionary stages of red giants in clus- 
ters of different ages. The depletion of the cluster sequence after 
logt w 8.5 can be explained by either cluster mortality (the typ- 
ical li fetime of local cluster s is only around logt ks 8.5, accord- 
ing to iPiskunov et al . 2006) or selection effects, which are more 
significant for older clusters. 

In the right column of Fig. [TJ we compare the distributions 
of the observed cluster colours with th e GALEV predictions. W e 
took the ages of the real clusters from lKharchenko et al.l (|2009), 
and calculated the integrated colours from the GALEV model. 
The discrepancy between the observed colours and those from 
the standard model is obvious. The local star clusters in our 
Galaxy are much bluer (by about 0.3 mag in the optical and 
about 0.8 mag in the near IR) than "theoretical" clusters of the 
same age. We discuss this below in connection with the results 
for the realistic model. 



3 http://stev. oapd. inaf . it/cgi-bin/cmd 
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Figure 2. Observed colours of 650 Galactic open clusters com- 
pared to theoretical colours computed from realistic and stan- 
dard SSP models. The upper row is for B—V, the bottom one 
for J— K s . The left column shows "integrated colour vs. log£ " 
diagrams, the right one the distributions of colour indices. Dots 
show Galactic open clusters from our sample. Black dots are 
young clusters used for constructing the colour distributions, 
grey dots the older clusters. The curves are the model tracks 
constructed in this paper. The dashed curve is for the stan- 
dard model, crosses show the realistic one. The observed colour 
distributions are shown as filled cyan histograms, the realistic 
model distribution as the hatched red ones. 



3. The discreteness of the IMF and the SSP 
models 

As we have already mentioned, a possible explanation of the 
observed disagreement is the continuity of the IMF adopted in 
standard SSP models, which means that every small mass inter- 
val of the given isochrone emits light according to its T e f f and 
log g. The total flux that a "theoretical" cluster emits is a result 
of the integration over the isochrone, weighted by the given IMF. 
It is immediately obvious that luminosity then scales with mass, 
once the age (or the isochrone) is fixed. In the discrete approach, 
however, only points on the isochrone where stars of a respective 
mass are actually present contribute. The random realisations of 
the IMF have statistical fluctuations even for a smooth proba- 
bility density distribution (i.e., smooth IMF). At the high-mass 
end of the mass range, these fluctuations are particularly large 
because of low number statistics, and produce gaps in the mass 
spectrum. In the case of a single supergiant present in the cluster 
at a given instant of time - a case frequently occurring in our lo- 
cal Galactic clusters - the colour fluctuations might be especially 
strong. If the main sequence lifetime of the star closest in mass 
to the supergiant exceeds the supergiant's lifetime, then after the 
supergiant's death the less-massive star remains on the main se- 
quence, which makes the cluster considerably bluer than it was 
before. The sparser the stellar population of a cluster is (i.e., the 
lower its total mass) the stronger the colour fluctuations are in 
the course of its evolution. After the most massive stars die, the 
red giant population of a cluster reaches an equilibrium and the 
colour fluctuations of a cluster begin to weaken. 

Fig. [2] is the counterpart of Fig. [T] for a realistic model 
of the scenario described above. The solid line with crosses 



and the dashed line shown here correspond to the discrete-IMF 
and to the continuous-IMF regimes, respectively, of the same 
model. The total initial mass of the model is assumed to be 
M c = 1000 m 0, and mi — 0.08 m Q . These parameters were 
chosen to provide a fit to the colour distribution in the sim- 
plest case of our local young cluster population. The ultimate 
upper mass limit m u = 100 m Q is defined by the Padova grid 
of isochrones, and m u depends on both of the adopted M c and 
mi. For example, for M c — 1000 777,0 and mi = 0.01, 0.1, 
and 1 TO0, the upper mass limit to the simulated mass spectrum 
equals m u « 60, 90, and 100 mQ. In the latter case, m u is lim- 
ited by the adopted evolutionary grid. 

For the standard models, this choice of parameters is how- 
ever not important. The derived evolutionary track depends on 
neither M c nor m; . The upper mass limit affects only the early 
stages of the cluster evolution. The only parameter influencing 
the standard model is the IMF slope a, although this influence 
is not strong. To change substantially the track, one would need 
to increase the exponent a to 4.35. A flattening of the IMF with 
respect to the Salpeter value of the slope does not have a strong 
effect on the colour-age relation. We note that our continuous- 
IMF model agrees with the predictions of the GALEV code 
and, hence, produces colours redder than the observed colours 
at log/; < 8.5. 

The realistic model with the above parameters produces 
a mass spectrum consisting of 3532 stars. In the mass range 
5 m© > m > 0.3 TO©, however, the number of stars is only 
599. This is of the same size as the number of Plei ades members 
observ ed in this mass range, which according to iBelikov et al.l 
(1998) equals 780. Thus, the model appears to be comparable 
to the Pleiades. The number of model stars with m > 5 m Q is 
equal to 20. So, before the model cluster reaches the age of the 
Pleiades (\ogt w 8.1), it has produced 20 red supergiants. The 
twenty corresponding "RG-events" are clearly seen in Fig. [2] 
At all other times (until logi w 8.1), the cluster is seen as a 
"main-sequence cluster" and resides in the MS-domain of the 
diagram. If we increase the lower bound mi, the number of mas- 
sive stars increases and the fraction of "RG-events" slowly in- 
creases, especially at older ages. But even at mi = 1 777-0, the 
"MS-clusters" are present in the track up to \ogt w 8.2. On 
the other hand, decreasing mi diminishes the number of "RG- 
events". In principle, from the frequency of red outliers in the 
"colour versus logt" diagram, one can estimate indirectly the 
parameter 777,;. 

Although the above example provides only an explanation 
for the local cluster population of typical mass of order of 
10 3 7770, it nevertheless illustrates the effect of the discreteness 
on extragalactic clusters where the bulk of cluster data comes 
from analysing t he integrated light. For clus ters of mass of the 
order of 10 5 m Q Dolph in & Kennicuttl (120021) found that this ef- 
fect is rat her small. This findi ng was interpreted by some work- 
ers (e.g., [Bastian et al. 2005) as a jus tification for negle cting 
discreteness. Among other arguments, Bastian et al. (20051) re- 
fer to the large masses of the observed extragalactic clusters, 
and assume that the effect only produces a symmetric spread 
around the "true" sequence. However, observations also find 
extragalactic clusters with masses typ ical of Galactic clusters 
logM c = 2 - 5 (see iMora et~ai1l2009l) . Since Fig.|2]shows that 
the effect is highly asymmetric with respect to the standard SSP 
sequence, a study of the effect over the entire range of cluster 
masses is appropriate. 

To estimate the dependence of the discreteness effect on 
cluster age and mass, we have performed a series of Monte Carlo 
simulations. We considered a set of cluster models with masses 



spread over the range logA/ c = 2.0 — 6.5. For every model of 
a particular mass at the starting point, the specific realisation of 
the initial mass spectrum was sampled at random from the IMF. 
The evolution of the synthetic cluster was then followed up to 
logi < 8.3. At every time step, we compared the colours of the 
realistic model with its standard counterpart and computed the 
colour differences 5 = I(B—V) r — I(B—V) S . Here the sub- 
script r represents "realistic", and s "standard". For every M c , 
we compiled N — 50 evolutionary sequences of cluster models 
by considering different statistical realisations of the initial mass 
spectrum, and evaluated at every time step the corresponding av- 
erage difference in colour 

N 
i=l 

and its statistical uncertainty e computed to be 

N 

e 2 = ^-A) 2 /(iV-l)/iV. 

i=l 

In Fig. [3] we summarize the colour bias produced by neglecting 
the IMF-discreteness effect, in the dependence on both the age 
and mass of a star cluster. The left panel shows the dependence 
of A on cluster age. The curves are for different model cluster 
masses. The bars indicate the uncertainty e. The first conclusion 
from Fig.|3]is that there is an overall negative bias. As expected, 
the maximum effect is observed at log t » 7, when the first red 
stars appear in the upper, sparsely populated CMD of a cluster. 
At younger ages, the effect diminishes since both realistic and 
standard models are populated with main-sequence stars only, 
which have approximately the same colour indices at the up- 
per MS. At \ogt > 7, the effect slowly weakens, reflecting the 
CMD geometry (distance between the MS and Hayashi limit). At 
log t ks 8, | A| becomes lower than 0. 1 for all considered masses. 
With increasing cluster mass, the effect in general diminishes. 
It is interesting to note that even at logA/ c = 6.5 the effect is 
small, but significant, and at \ogt = 6.7 — 7.0 the colour shift 
is of the order of a few hundredths of a magnitude. For masses 
log M c > 4, |A| remains below 0.1 at all ages considered. 

The explicit dependence of the discreteness effect on cluster 
mass, providing the net-estimate of the effect for young clusters, 
can be derived if one averages the A versus log t relations over 
the entire range of ages. The result is shown in the right panel of 
Fig. |3] The effect is considerable (|(A) t | > 0.05) for logM c < 
4. At higher masses it is small although significant. 

For the other colours the behaviour is qualitatively similar, 
although the redder the passband, the stronger the effect is. This 
is intuitively understandable. The effect is caused by quantiza- 
tion, because red stars radiate proportionally more of their light 
at longer wavelengths. For example, in the case of I(J—K S ) the 
largest difference between the realistic and the standard models 
is about —0.9 (cf. Fig. [3j. In a forthcoming paper, this issue will 
be discussed in more detail. 

4. Concluding remarks 

Observations of local Galactic open clusters have enabled us to 
measure their integrated magnitudes, masses, ages, and redden- 
ing. Using this data set, we have found a remarkable discrepancy 
between the observed colours and the predictions of SSP mod- 
els. The main reason for this disagreement is the neglect of the 
assumption of IMF-discreteness. When this effect is taken into 
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Figure 3. The colour displacement versus age relations (left) for 
cluster models of different masses. The curves correspond (from 
bottom to top) to logM c = 2, 3, 4, 6.5. The vertical bars indi- 
cate the statistical uncertainty e. In the right panel, we show the 
displacement averaged over the age interval logt = 6.5 — 8.3 
as a function of cluster mass. The bars indicate the error in the 
averaging. 

account, the model agrees adequately with the observations and 
is even able to explain the large colour spread observed in the 
empirical colour-age relation in a natural way. 

In which conditions is the effect of discreteness relevant? 
Since it is a consequence of low-number statistics, it depends 
on the sparseness of the stellar population in the upper CMD of 
a cluster where the bulk of light is emitted. The density of the 
population depends primarily on cluster mass, and in addition 
on cluster age, on [strange at first glance] the lower mass limit 
of the stars formed, and to some degree on the slope of the IMF. 
In the context of this letter, other factors are not important. 

According to our investigation, the systematic offset between 
the continuous- and discrete-IMF colours diminishes substan- 
tially but not completely at \ogt > 7.9, at cluster masses M c > 
1O 5 TO0. At younger ages, it remains present even in massive 
clusters, and for M c ^ 10 4 it is larger than 0.1 mag in (B—V). 
Only for very massive clusters (M c > 10 6 m@) and young ages 
(logf < 7.5), the offset falls below a typical observational error. 
We note in passing that the effect is stronger for redder pass- 
bands. The se findings are in good agree ment with the theoretical 
forecast of ICervifio & Luridia na (2004). The immediate conse- 
quence of the application of a continuous-IMF approach to the 
SSP models of stellar clusters is a systematic underestimate of 
reddening. Other problems affect the age and mass determina- 
tion based on evolutionary grids of these models (at least for 
masses below 10 6 m Q ). Because of the rather flat colour-age re- 
lation of young clusters, a systematic error of the order of 0.1 
mag produces an error in cluster age of about one order of mag- 
nitude. On the other hand, there is a danger that the existence of 
a large number of "main-sequence" (i.e., blue) clusters could be 
interpreted as evidence of a recent burst of star formation. We 
can exclude such a burst having occured in our Galactic open 
clusters even if they are blue. Finally, the technique of the param- 
eter determination should be changed to incorporate the cluster 
flux variations caused by the discrete nature of the upper IMF. 
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